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We estimate the number of templates, computational power, and storage required for a one-step 
matched filtering search for gravitational waves from inspiraling compact binaries. Our estimates 
for the one-step search strategy should serve as benchmarks for the evaluation of more sophisti¬ 
cated strategies such as hierarchical searches. We use a discrete family of two-parameter waveform 
templates based on the second post-Newtonian approximation for binaries composed of nonspinning 
compact bodies in circular orbits. We present estimates for all of the large- and mid-scale inter¬ 
ferometers now under construction: LIGO (three conhgurations), VIRGO, GEO600, and TAMA. 
To search for binaries with components more massive than mmin = 0.2Mq while losing no more 
than 10% of events due to coarseness of template spacing, the initial LIGO interferometers will 
require about 1.0 x 10^^ flops (floating point operations per second) for data analysis to keep up 
with data acquisition. This is several times higher than estimated in previous work by Owen, in 
part because of the improved family of templates and in part because we use more realistic (higher) 
sampling rates. Enhanced LIGO, GEO600, and TAMA will require computational power similar 
to initial LIGO. Advanced LIGO will require 7.8 x 10^^ flops, and VIRGO will require 4.8 x 10^^ 
flops to take full advantage of its broad target noise spectrum. If the templates are stored rather 
than generated as needed, storage requirements range from 1.5 x 10^^ real numbers for TAMA to 
6.2 X 10^'* for VIRGO. The computational power required scales roughly as and the storage 

as Since these scalings are perturbed by the curvature of the parameter space at second 

post-Newtonian order, we also provide estimates for a search with mmin = IMq. Finally, we sketch 
and discuss an algorithm for placing the templates in the parameter space. 

PACS numbers; 04.80.Nn, 07.05.Kf, 97.80.-d 


I. INTRODUCTION 

Close binary systems composed of compact objects (such as black holes and neutron stars) are expected to be 
an important source of gravitational waves for broadband laser interferometers such as LIGO, VIRGO, GEO600, 
and TAMA The orbit of a compact binary decays under the influence of gravitational radiation reaction, 

emitting a gravitational wave signal that increases in amplitude and “chirps” upward in frequency as the objects 
spiral in toward each other just before their final coalescence. According to current astronomical lore i-0 , the rate 
of coalescences should be about three per year within 200 to 300 Mpc Q of the Earth for neutron-star-neutron-star 
binaries, and within 400 Mpc to 1 Gpc for black-hole-black-hole binaries. Signals from inspiraling compact binaries 
at these distances are strong enough to be detected by “enhanced” LIGO ||^ and VIRGO interferometers, but only if 
aided by a nearly optimal signal-processing technique. Fortunately, the distinctive frequency chirp has been calculated 
to a remarkable degree of precision using a variety of approximations to the general relativistic two-body problem 
(e.g. P,|T^). Because the functional form of the chirp is quite well-known, a search for inspiral signals in noisy data 
is ideally suited to matched filtering. 

Matched hltering njlhas long been known to be the optimal linear signal-processing technique and is well-discussed 
in the literature (e.g. ||l2| ); therefore we will only briefly summarize it here. In the frequency domain, a matched hlter 
is a best-guess template of the expected signal waveform divided by the interferometer’s spectral noise density. The 
interferometer output is cross-correlated with the matched hlter at different time delays to produce a hltered output. 
The signal-to-noise ratio, dehned as the ratio of the actual value of the hltered output to its rms value in the presence 
of pure noise, is compared to a predetermined threshold to decide if a signal is present in the noise. If the signal from 
which the matched hlter was constructed is present, it contributes coherently to the cross-correlation, while the noise 
contributes incoherently and thus is reduced relative to the signal. Also, the weighting of the cross-correlation by 
the inverse of the spectral noise density emphasizes those frequencies to which the interferometer is most sensitive. 
Consequently, signals thousands of cycles long whose unhltered amplitude is only a few percent of the rms noise can 
be detected. 
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A matched filtering search for inspiraling compact binaries can be computationally intensive due to the variety of 
possible waveforms. Although the inspiral signals are all expected to have the same functional form, this form depends 
on several parameters—the masses of the two objects, their spins, the eccentricity of their orbit, etc.—some of them 
quite strongly. A filter constructed from a waveform template with any given parameter vector may do a very poor 
job of detecting a signal with another parameter vector. That is, the difference in parameter vectors can lead to a 
greatly reduced cross-correlation between the two wave forms; and in general, the greater the difference, the more the 
cross-correlation is reduced. Because the parameter vector of a signal is not known in advance, it is necessary to filter 
the data with a family of templates located at various points in parameter space—e.g., placed on a lattice—such that 
any signal will lie close enough to at least one of the templates to have a good cross-correlation with that template. 

There are several questions that must be answered in order to determine the feasibility of a matched filtering search 
strategy and, if feasible, to implement it. Which parameters significantly affect the waveform? How should the spacing 
of the template parameters (lattice points) be chosen? Is there a parametrization that is in some sense “preferred” by 
the template waveforms? How many templates are needed to cover a given region of interest in the parameter space, 
and how much computing power and memory will it cost to process the data through them? In the case of inspiraling 
compact binaries, the full general-relativistic waveforms are not exactly known, but are instead approximated (e.g., 
using the post-Newtonian scheme); and we must also ask, what approximation to the true waveforms is good enough? 

All of these questions have been addressed in recent years, at least at some level. The current state of affairs is 
summarized by the following brief review of the literature: 

The standard measure for deciding what class of waveforms is good enough is the fitting factor {FF) introduced by 
Apostolatos l]^. The fitting factor is effectively the fraction of optimal signal-to-noise-ratio obtained when filtering 
the data with an approximate family of templates. Because binaries are (on large scales) uniformly distributed in space 
and because the signal strength scales inversely with distance, the fraction of event rate retained is approximately 
FF^. Therefore it has become conventional to regard FF = 97%—i.e., 10% loss of event rate—as a reasonable goal. 
Using the standard post-Newtonian expansion in the test-mass case (i.e., when one body is much less massive than 
the other so that the waveforms can be computed with arbitrarily high precision using the Teukolsky formalism), Droz 
and Poisson 0 found that 2PN signals had fitting factors of 90% or higher. Damour, Iyer, and Sathyaprakash 0] 
have devised a new way of rearranging the usual post-Newtonian expansion (similar to the way Pade approximants 
rearrange the coefficients of a Taylor expansion) to take advantage of physical intuition in constructing templates. 
They find fitting factors of 95% or higher for the 2PN templates. Research underway by the authors of will lead 
to 3PN templates that should easily achieve FF > 97%. 

Several people |^,|^,|^ have shown that it is insufficient to use templates that depend on just one shape parameter 
(the “chirp mass,” which governs the rate of frequency sweep at Newtonian-quadrupole order). To achieve FF > 90% 
one must include the masses of both objects as template parameters, as was done in the above 2PN analyses [ p^ , p^ , 
and as is being done in the forthcoming 3PN templates. 

Apostolatos showed that, for binaries whose components spin rapidly about their own axes which are 

orthogonal to the orbital plane so that there is no precession, neglecting the spin parameters (i.e., using two-mass- 
parameter waveforms based on the theory of spinless binaries) degraded the fitting factors by less than 2%. With 
precession the situation is much more complicated, and data analysis algorithms are as yet poorly developed: It is clear 
that there are interesting corners of parameter space (most especially a neutron star in a substantially nonequatorial, 
processing orbit around a much more massive, rapidly spinning black hole) in which the two-mass-parameter spinless 
waveforms give FF 90%; to search for such binaries will require waveforms with three or more parameters [ ^ . 
However, the 2PN (or 3PN) two-mass-parameter waveforms do appear to cover adequately a significant portion of 
the parameter space for processing binaries [|^ . 

Sathyaprakash showed that in computations with the two-mass-parameter waveforms, the best coordinates to 
use on the parameter space are not the two masses, but rather the inspiral times from some fiducial frequency to final 
merger, as computed at Newtonian and first post-Newtonian order. Working with the restricted first post-Newtonian 
wave forms (see below) he found that the effective dimension of the parameter space is nearly one. 

Sathyaprakash and Dhurandhar ||20|-|^ developed a criterion for putting templates at discrete points on a grid in 
parameter space and numerically implemented their criterion for a one- par ameter (Newtonian) family of templates and 
for simple noise models. They introduced the concept of what Owen [|23|| later called the minimal match (analagous 
to the fitting factor) as a measure of how well a discrete set of templates covers the parameter space and estimated 
the computational costs for an on-line search. 

Owen [ p^ , building on the work of Sathyaprakash and Dhurandhar, defined a metric on the parameter space from 
which one can semi-analytically calculate (i) the template spacing needed to achieve a desired minimal match, (ii) the 
total number of templates needed, and (iii) the computational requirements to keep up with the data—for any family 
of waveforms and any interferometer noise spectrum. Owen combined this metric-based formalism with computational 
counting procedures from Schutz to estimate the computational requirements for LIGO searches based on two- 
parameter IPN templates. These estimates were confirmed by Apostolatos ||^ using a numerical method in the 
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vein of (but more sophisticated than) the previous work of Sathyaprakash and Durandhar ||20|-|22| . Apostolatos also 
showed that a search for processing binaries that fully covers all the nooks and crannies of the precessional parameter 
space, using currently available templates and techniques, is prohibitively costly. 

Mohanty and Dhurandhar |2|,^ have studied hierarchical search strategies. Such strategies reduce computational 
costs by making a first pass of the data through a coarsely-spaced template grid and a low signal-to-noise threshold 
to identify candidate signals. Each candidate flagged by the first pass is examined more closely with a second, 
finely-spaced grid of templates and a higher threshold to weed out false alarms. Such strategies can reduce the total 
computational requirements by roughly a factor 25. 

The purpose of this paper is to refine and update the analyses by Owen for the two-parameter, spinless 
templates that are likely to be used for binary-inspiral searches in ground-based interferometers. This refinement is 
needed because the kilometer-scale interferometers will begin taking data in about 2 years (preliminary, engineering 
run); people are even now designing software to implement the simplest matched filtering search algorithm; and in the 
context of these implementations, the factor of 3 accuracy attempted in Ref. is no longer adequate. The numbers 
that are derived in this paper should establish a reliable baseline cost to which more sophisticated search strategies 
(e.g., hierarchical searches) can be compared. 

The substantial differences between this paper and Ref. are that we now (i) approximate the phase evolution 
of the inspiral waveform to 2PN rather than IPN order; (ii) give results for the noise spectra of several more inter¬ 
ferometers; and (iii) use a better estimate of the sampling frequency needed for each interferometer. We assume the 
following fiducial search: a minimal match of 97% (corresponding to 10% loss of event rate due to coarse parameter 
space coverage), second post-Newtonian waveforms, and templates made for objects of minimum mass rumin = Q.2Mq 
and up. 

Our results for the computational requirements are given in Tables These tables show that the initial LIGO 

interferometers need about twice as many templates and triple the computational power estimated in Ref. . These 
increases result mainly from using 2PN waveforms rather than the (clearly inadequate) IPN, and from a using a 
higher sampling rate (as, it turns out, is required to keep time-step discretization error from compromising the 97% 
minimal match). GEO600 requires slightly more templates and power than LIGO because of its flatter noise spectrum, 
while TAMA requires slightly less because its sensitivity is limited to higher frequencies where there are fewer cycles. 
Initial VIRGO, with its extremely broad and flat spectrum, requires about the same as advanced LIGO. We conclude 
that even if VIRGO achieves its targeted low-frequency sensitivity, it will be unable to take full advantage of that 
sensitivity with matched filtering because of the large number of templates required. 

The rest of this paper is organized as follows. In Sec. II we analyze the application of matched filtering to a 
search for inspiraling binaries and summarize the method of Ref. |23| which uses differential geometry to answer 
important questions about such a search. We use this method in Sec. Ill to estimate the computational costs and 
other requirements of a matched-filtering binary search for LIGO, VIRGO, GEO600, and TAMA. In Sec. IV we 
illustrate a detailed example of a template placement algorithm, and in Sec. V we discuss our results. 


II. FORMALISM 


This Section summarizes material previously presented in with several incremental improvements. We begin 
by introducing some notation. 

The Fourier transform of a function h{t) is denoted by h{f), where 

/ OO 

( 2 . 1 ) 

-OO 


We write the interferometer output h{t) as the sum of noise n{t) and a signal As{t), where we have separated the 
signal into a dimensionless, time-independent amplitude A and a “shape” function s{t) which is defined to have unit 
norm [see Eq. (2.4) below]. 

The strain power spectral noise density of an interferometer is denoted by Sh{f)- We use the one-sided spectral 
density, defined by 


E[h(/i)r(/2)] = - /2)5,(|/i|), (2.2) 

where E[ ] denotes the expectation value over an ensemble of realizations of the noise and an asterisk denotes complex 
conjugation. 

We use geometrized units, i.e. Newton’s gravitational constant G and the speed of light c have values of unity. 
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A. Matched filtering 


First we flesh out the Introduction’s brief description of matched filtering. In the simplest idealization of matched 
filtering, the filtered signal-to-noise ratio is defined by 


S _ {h,u) 

N rms (n, u) 


(2.3) 


Here u is the template waveform used to filter the data stream h, and the inner product 


(a, b) = 4 Re 



Shif) 


(2.4) 


is the noise-weighted cross-correlation between a and b (cf. |^). The denominator of (|2.3|) is equal to {u, w), the 
norm of u (see Sec. II B of Ref. for a proof). Because the norm of u cancels out of (|2.3|) , we can simplify our 
calculations without loss of generality by considering all templates to have unit norm. 

When searching for a parametrized family of signals the situation is somewhat more complicated. The parameter 
values of the signals are not known in advance; therefore one must filter the data through many templates constructed 
at different points in the parameter space. To develop a strategy for searching the parameter space, one must know 
how much the S/N is reduced by using a template whose parameter values differ from those of the signal. Neglecting 
fluctuations due to the noise, the fraction of the optimal S/N obtained by using the wrong parameter values is given 
by the ambiguity function, 


A{X,A) = {u{X),u{A)) (2.5) 

(see, e.g., Chs. XIII and X of Ref. @)- Here A and A are the parameter vectors of the signal and template (it 
doesn’t matter which is which). The ambiguity function, as its name implies, is a measure of how distinguishable two 
waveforms are with respect to the matched filtering process. It can be regarded as an inner product on the waveform 
parameter space and is fundamental to the theory of parameter estimation |T^J^ . 

For the purposes of a search for inspiraling compact binaries, the ambiguity function isn ’t q uite what is needed. 
This is because the test statistic (for a given set of parameter values 6) is not given by Eq. (|2.3| ), but rather by 

max ^- - -—-— 2.6 

(j}a,tc rms (n, u[U)) 

Here (j>c and t^ are respectively the coalescence time and coalescence phase. We separate these parameters out from 
the rest: A = {(j}i,,tc,d), where 9 is the vector of intrinsic parameters that determine the shape of the waveform and 
4>c and tc are extrinsic parameters ||^ (also refereed to as kinematical and dynamical parameters respectively). 

The practical difference is that maximization over the extrinsic parameters is perf ormed automatically by Fourier 
transforming, taking the absolute value, and looking for peaks. The use of Eq. ( ^. 6 [ ) as a detection statistic suggests 
the definition of a modified ambiguity function known as the match 

M(ei, 6 / 2 ) = max (u( 6 >i), u( 6 > 2 )e*( 2 -/*=-^=)), (2.7) 

where the templates u are assumed to have unit norm. The use of this match function rather than the ambiguity 
function takes into account the fact that a search can benefit from systematic errors in the extrinsic parameters. 


B. Applications of differential geometry 


The match (2.7) can be regarded as an inner product on the space of template shapes and intrinsic template 


parameters, and correspondingly one can define a metric on this space 


^ _i d^M{e,&) 

’ 2 de^dQ3 


Qk^gk 


( 2 . 8 ) 


The metric (2.8) is derived from the match © in the same way the information matrix 


ambiguity function 


is derived from the 

and plays a role in signal detection similar to that played by the information matrix in 
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parameter estimation. The gij can be derived by expanding M{6, ©) about ® = 6, or equivalently by projecting r^ 
on the subspace orthogonal to (f>c and tc- 

The gij can be used to approximate the match in the regime 1 — M <C 1 by 


M(e, e + Ae) ~ 1 - A6»*A6»^ 


(2.9) 


which is simply another way of writing the Taylor expansion of M{9, 9 + A6) about A9 = 0. (The first derivativ e te rm 
vanishes because M takes its maximum value of unity at A9 = 0.) We find that the quadratic approximation ( ^.9| ) is 
good typically for M ~ 0.95 or greater, though this depends on the waveform and noise spectrum used. Experience 
suggests that the quadratic approximation gene rally underestimates the true match; and thus the spacings and 
numbers of templates we calculate using Eq. (|2.9| ) err on the safe side. See Fig. for an example. 

In the limit of close template spacing, Eq. ( |2.9D leads to a simple, analytical way of placing templates on a lattice. 
We discuss this in some detail in Sec. fV, but for now turn our attention to the use of the quadratic approximation 
in calculating the number of templates needed for a lattice. 


C. Computational costs 


If the number M of templates needed to cover a region of interest is large, it is well approximated by the ratio of 
the proper volume of the region of interest to the proper volume per template AE, 


JV = {AV)-^ J d^0ijdet\\g.j\\, 


( 2 . 10 ) 


where D is the dimension of the parameter space |^. Equation ( 2.10 ) underestimates Af when not in the limit 
AV ^ 0 ( Af —> oo). The reason is template spill over, i.e. the fact that in any real algorithm for laying out 
templates, those on the boundaries of the region of interest will to some extent cover regions just outside. This effect 
is small in the limit of many templates because it goes as the surface-to-volume ratio of the region of interest. 

The proper volume per template, AV, depends on the packing algorithm used, which in turn depends on the 
number D of dimensions (see Sec. |^. For D = 2, the optimal packing is a hexagonal lattice, and thus 


AE= —(1-MM), 


( 2 . 11 ) 


where MM is the minimal match parameter defined in Ref. [ P3[ as the match between signal and template in the case 
when the signal lies equidistant between all the nearest templates (i.e., the worst-case scenario). There is no packing 
scheme which is optimal for all D, but it is always possible (though inefficient) to use a hypercubic lattice, for which 


AE = {2iJ{l- MM)/D)°. 


( 2 . 12 ) 


For inspiraling compact binaries. Ref. |P3| ha s spelled out a detailed prescription for obtaining the needed to 


evaluate the proper volume integral in Eq. (2.10|). In summary, first one obtains a metric including the tc parameter, 


lap = ■;r(J"[V'aV'/3] - J[^a]J[^p\), 


(2.13) 


where ipct is the gradient of the waveform phase dt in the parameter space of intrinsic parameters plus tc and the 
moment functionals 


J[a] = 


(r^/^a(/)) 

4/o“rf/^ 


(2.14) 


JO ShU) 

can be expanded (for binary chirp waveforms) in terms of the noise moments ]29 


J{p) = 


(///o)_^ 

L.o ^ Sn{f) 


-1 


/o Snif) 


(2.15) 


where /o is the frequency of the minimum of Sh{f)- Then one projects out the coalescence time tc to obtain 

gij = lij - loilojhm- (2.16) 
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FIG. 1. Comparison of the match to the quadratic approximation in the case of: (a) first post-Newtonian waveforms and 
(b) second post-Newtonian wave forms. The solid line is obtained using the analytical model Eq. (2.9) and the dots denote the 
match obtained by numerically computing the ‘distance’ from the centre, at which the match drops to the chosen value 0.97. 
The two agree very well in the case of first post-Newtonian waveforms while in the second post-Ne wton ian case the agreement 
is quite good for all but waveforms of massive binaries. The coordinates (ro,ri) are defined in Eq. (3.2). 


6 










Once Af has been found it is a relatively straightforward matter to calculate the CPU power and storage required 
to process all the templates in an on-line search. The interferometer data stream is broken up for processing into 
segments of D samples (real numbers), such that D ^ F where F is the length (in real numbers) of the longest 
filter. (See Schutz for a discussion of the optimization of D/F, taking into account the fact that successive data 
segments must overlap by at least F to avoid circular correlations in the Fourier transform.) Using the operations 
count for a real Fourier transform |^] , filtering the data segment through Af templates of length F requires 

Af D{W + 3\og2F) (2.17) 

floating point operations. If we take the sampling frequency to be 2/^ (see Sec. Ill and Table |), the computational 
power required to keep pace with data acquisition is 

IP ~Af/432 + 61og2F) (2.18) 

flops (floating point operations per second). 


III. COMPUTATIONAL COST USING RESTRICTED 2PN TEMPLATES 

In this section, using the geometric formalism summarized in Sec. II, we calculate the number A/" of templates 
required to cover a region of interest as a function of the minimal match. We then use this number to calculate the 
computational cost of filtering a single interferometer’s output through all these templates in an on-line search. 


A. Functional form of the templates 


We construct our waveform templates using two intrinsic parameters based on the masses of the binary’s com¬ 
ponents. Inspiral waveforms in principle can be strongly affected by several other parameters: spins of the two 
components, orbital eccentricity, and several angles describing the orientation of the binary with respect to the inter¬ 
ferometer. However, it is believed that two-parameter templates will be adequate to search for most binaries for the 
following reasons. 

(i) Based on models of the evolution of currently known binary pulsars, it is expected that typical NS-NS binaries 
will have spins of negligible magnitude (spin/mass^ ^ 1). Apostolatos has shown that, even if the magnitudes of 
the spins are large, their effect on S/N is small (reduces FF by less than 2% if the orbit and the spin vectors do not 
process). He has also shown that precession will not reduce FF below 90% except in the relatively small region of 
parameter space containing binaries with a neutron star orbiting a more massive, rapidly rotating (spin/mass^ ~ 1) 
black hole with orbital angular momentum inclined by more than about 30 degrees to the black hole’s spin. 

(ii) It has long been known that gravitational radiation reaction circularizes all but the most eccentric orbits 
on a time scale much smaller than the lifetime of the binary if the progenitor system was the same binary. This may 
not be true, however, in the case of close binaries formed by capture in densely populated environments, e.g. galactic 
nuclei/globular clusters. 

(iii) The angles make no difference in our analysis because we use the restricted post-Newtonian approximation , 
in which the phase evolution of the inspiral waveform is followed to a high post-Newtonian order but the amplitude 
is only followed to lowest order. In this approximation, the combined effect of the angles is to multiply the waveform 
by a constant amplitude and phase factor, which does not affect the choice of search templates Presently it is 
believed that the restricted post-Newtonian approximation will be good enough for data analysis of ground-based 
interferometers. 

The standard post-Newtonian expansion of the waveform phase is given as a function of mass parameters based on 
the standard astronomical choices M (total mass) and p, (reduced mass). In order to clearly isolate test-mass terms 
(i.e., those that remain when one body is much less massive than the other), the symmetric mass ratio rj = p/M is 
typically used instead of p. In terms of M and 77 , the second post-Newtonian waveform phase can be calculated from 
the energy loss formula of Blanchet et al. |p| as 


-q) 


128 


{irMf) 


- 5/3 


-flO 


3 058 673 
1016 064 



5 429 
1008 


- 167r(7rM/) 


(3.1) 
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[cf. Eq. (3.6) of Poisson and Will |p^ ]. However, M and r] are inconvenient parameters for our purposes because, when 
they are used as the parameter-space coordinates, the values of the metric components vary strongly over parameter 
space, making calculations unnecessarily difficult and prone to numerical error. 

In earlier analyses it was found more convenient to use as parameters the Newtonian and IPN chirp times 


ro = 

5 , ,_2 ( 743 11 

'^1 = (^/o) + 


192 


V3367? 


(3.2a) 

(3.2b) 


which are respectively the Newtonian and IPN contributions to the time it takes the instantaneous gravitational- 
wave frequency to (formally) evolve from /o to infinity. The chirp times are more convenient than the usual mass 
parameters because, when they are chosen as parameter-space coordinates, at IPN order the metric components 
are constant. Assuming the post-Newtonian expansion has reasonable convergence properties, one would expect the 
metric components in these coordinates to remain nearly constant at higher post-Newtonian orders (and indeed we 


find this is so). However, at higher than IPN order one cannot write the waveform phase analytically in terms of tq 
and Ti. To remedy this, following Mohanty |^|, we base our second parameter on the 1.5PN chirp time ri.s (see 
for a definition). More specifically, we introduce new dimensionless coordinates in parameter space. We define 


= 27r/oro = \ 

= 27r/ori,5 = ^(7rM/o)“^/3^-i^ 

which can be inverted to obtain 


]\/r 

5 

327r2/o 9^ ’ 


■I67r3 (6»i) 

7 ] — 

25 (02) 


(3.3a) 

(3.3b) 


(3.4a) 

(3.4b) 


This choice of (0^, 9^) lets us write the waveform phase analytically while keeping the metric components from varying 
too strongly; therefore it is convenient for calculating numbers of templates (see below). However, for other purposes 
(tq, Ti) are just as convenient, and to be consistent with the literature we will use them. 


B. Noise spectra 

In this paper we consider the noise spectra of the four large- and intermediate-scale interferometer projects, LIGO, 
VIRGO, GEO600, and TAMA. For LIGO we use three noise spectra corresponding to three interferometer configu¬ 
rations, the “first interferometers” Q (which are planned to perform a gravitational-wave search in 2002-2003), the 
“enhanced interferometers” (which are likely to be carrying out searches in the mid 2000’s), and the “advanced 
interferometers” Q (which are thought representative of the type of detector that might operate a few years after 
the enhanced ones). For convenience we abbreviate these, respectively, as LIGO I, LIGO H, and LIGO HI. We use 
the VIRGO and GEO600 noise spectra given in Refs. |^,||. For TAMA we use the noise spectrum given by M.-K. 
Fujimoto in Ref. li¬ 
lt is desirable to have simple analytical fits to the noise power spectral densities used. Kip Thorne and Scott Hughes 
(private communications) have provided us with fits to LIGO I and LIGO HI, respectively, and a fit to LIGO H was 
derived in Ref. We have constructed our own fits to the remaining noise spectra listed in the previous paragraph. 
All of these analytical fits to the noise spectra are tabulated in Table ^ 

The shape of each noise spectrum determines natural low- and high-frequency cutoffs for the matched filtering 
integrals. The low-frequency cutoff fs is the defined as the frequency above which 99% of S'^/N^ is obtained; we 
call this the seismic frequency and denote it with a subscript s because it is typically near the frequency at which 
seismic noise causes the noise power spectral density Sh{f) to begin rising sharply. Because integration below this 
frequency contributes little to the detectability of signals and costs much in terms of total number of templates and 
computational resources, we assume the templates to begin with gravitational-wave frequency fs- 
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TABLE I. Analytical fits to noise power spectral densities Sh{f) of the interferometers treated in this paper. Here So is the 
minimum value of Sh{f), and fo is the frequency at which the minimum value occurs. For our purposes Sh{f) can be treated 
as infinite below the seismic frequency fa. The high-frequency cutoff /„ is chosen so that the loss of signal-to-noise ratio due 
to finite sampling rate 2fu is 0.75% (see text). 


Detector 

Fit to noise power spectral density 

So (Hz-1) 

fo (Hz) 

/. (Hz) 

fu (Hz) 

LIGO I 

So/S [ifo/f)^ + 2(///o)2] 

8.0 X IQ-i® 

175 

40 

1300 

LIGO II 

So/ll {2 (/o//)®/ 2 +9/2[l + (///o)=]} 

7.9 X 10-1® 

110 

25 

900 

LIGO III 

So/5 {(/o//)" + 2[l + (///o)^]} 

2.3 X lO-i® 

75 

12 

625 

VIRGO 

So/4: [290(/.//)5 + 2(/o//) + 1 + (///o)2] 

1.1 X 10-1® 

475 

16 

2750 

GEO600 

So/5 [4 (/o//)^/2-2 + 3(///o)2] 

6.6 X 10-1® 

210 

40 

1450 

TAMA 

So/32 {(/o//)® + 13(/o//) + 9[1 + if/fo)^]} 

2.4 X 10-11 

400 

75 

3400 


The high-frequency cutoff /„, as pointed out by Eanna Flanagan (private communication), needs to be high enough 
that the S/N degradation due to discrete time steps At = l/(2/„) in the data analysis is less than that due to discrete 
choices of the templates’ intrinsic parameters. Quantitatively, this requires that 


7oo/(4/„)' < 1 - MM, (3.5) 

where 700 is the t^-tc component of the metric before is projected out [see Ref. 1^, Eq. (2.23)]. At any post- 
Newtonian order we have 700 = (27r/o)^[J(l) — J(4)^]/2 [see Ref. [^, Eq. (2.27)] and thus 


fu > tt/o 


j(i) - 

8(1 - MM) ■ 


(3.6) 


We have chosen /„ to be twice the right-hand side of this expression so that the loss due to sampling is 1 /4 the loss 
due to discrete values of the intrinsic parameters. There is probably a more clever way of optimizing /„, but this is a 
first cut. 


C. Number of templates and computational cost 


We now proceed with the calculation of the n umb er of templates needed to perform a single-pass, on-line search 
for gravitational-wave signals of the form in Eq. (^), for the noise spectra in Table |. 

The first task is to obtain the intrinsic-parameter metric. We rewrite the waveform phase at 2PN order as 




- 5/3 
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(3.7) 


From Eqs. (13) and Eqs. (|2.13|) -( ^.16|) it is straightforward to derive the metric components gij with symbolic 
manipulation software. However, the general expressions for the gij (and even the 7 ij) are too lengthy and complicated 
to display here. 

Next, we obtain the proper volume 


R = J ¥e ^detll 5 ,,ll ( 3 . 8 ) 

of the region of interest. The boundaries of this region are set by the range of masses [mmin, Wmax] of the individual 
objects in binaries. It turns out that only rumin has a strong influence on V ]p0| , p^ , p^ . For IPN waveforms the intrinsic 
parameter space is flat (i.e. det \\gij jj is constant) and thus V can be obtained analytically using the coordinates (tq, ti). 
Beyond IPN order in the waveform, the intrinsic parameter space is not flat, and thus it is convenient to calculate 
the proper volume numerically, e.g. by a Monte Carlo method. In order to make the integrand as nearly constant as 
possible we use the metric tensor transformation law to switch from (0^,0^) coordinates to (ro,Ti) [because in these 
coordinates the metric components vary more slowly over the region of interest than in (0^,02)]. Finally, to obtain 
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TABLE II. Numbers of templates required to cover parameter space at a minimal match of 97% with a rectangular lattice. 
The region of interest is that inhabited by binaries with component masses greater than mmin. We use restricted post-Newtonian 
(PN) templates whose phase evolution is accurate to the indicated order. 


Detector 

IPN 

rrimin = 0.2Mq 

1.5PN 

2PN 

mmin = IM0 

2PN 

LIGO I 

2.3 X 10“ 

5.6 X 10® 

4.8 X 10“ 

1.1 X 10'* 

LIGO II 

1.0 X 10® 

1.9 X 10® 

1.7 X 10® 

4.0 X 10* 

LIGO III 

5.7 X 10® 

7.9 X 10® 

6.9 X 10® 

1.7 X 10® 

VIRGO 

5.8 X 10® 

1.1 X 10 '^ 

9.3 X 10® 

2.2 X 10® 

GEO600 

4.2 X 10® 

9.7 X 10® 

8.3 X 10® 

1.9 X 10* 

TAMA 

5.1 X 10'* 

1.7 X 10® 

1.4 X 10® 

3.1 X 10® 


TABLE III. Computational costs obtained from the numbers of templates given in Table ^ using 2PN templates and 
m-rain = Q.2Mq. The symbol F denotes the length (in real numbers) of the longest template. 


Detector 

logzF 

CPU power P (flops) 

Storage S (reals) 

LIGO I 

21 

9.9 X 10*® 

1.0 X lO*"" 

LIGO II 

22 

2.5 X 10** 

7.1 X 10*° 

LIGO III 

25 

7.8 X 10** 

2.3 X 10** 

VIRGO 

26 

4.8 X 10*^ 

6.2 X 10** 

GEO600 

21 

1.9 X 10** 

1.8 X 10*° 

TAMA 

20 

7.3 X 10*° 

1.5 X 10** 


the number of templates we divide the proper volume V b y the proper volume per template AI/. For a rectangular 
lattice in two dimensions, AI^ = 2(1 — MM) [cf. Eq. ( ^.12| )]. In Table || we give values of AC calculat ed fo r the noise 
spectra given in Table | assuming a rectangular lattice. For a hexagonal lattice, AE is given by Eq. ( ^.11 ), resulting 
in a reduction of about 20% in AC, V, and S. However, in practice much of this reduction may be offset by the 
details of actually constructing a lattice (see Sec. IV) and therefore we make our estimates for the more conservative 
rectangular case. 

These numbers are easily translated into the computational costs shown in Tables HI and IV: To a good approx¬ 
imation, the length of the longest filter (with mi = m 2 = mmin, i-e. equal mass binaries so that 77 = 1/4) is given 
by 


F ~ 2/„ro(/o//s)®/3 


-5/3 


The storage required for all of the templates is then roughly given by 

S-NF, 


(3.9) 


(3.10) 


and the computational power V is given by Eq. (2.18). 


D. Scaling laws 

Approximate scaling laws can be obtained as follows. At IPN order the intrinsic parameter space is flat, and thus 
the proper volume of the region of interest scales as the (dimensionless) coordinate volume, i.e. 


TABLE IV. 


Computational costs as in Table III 


except here we assume rrimin = IMq. 


Detector 

logaF 

CPU power P (flops) 

Storage S (reals) 

LIGO I 

17 

1.9 X 10° 

1.4 X 10° 

LIGO II 

19 

5.2 X 10° 

2.1 X 10*° 

LIGO III 

21 

1.7 X 10*° 

3.6 X 10** 

VIRGO 

22 

9.8 X 10*° 

9.1 X 10** 

GEO600 

17 

3.7 X 10° 

2.5 X 10° 

TAMA 

16 

1.3 X 10° 

2.0 X 10® 
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V ~ foTo X foTi ~ 


( 3 . 11 ) 


At higher post-Newtonian orders the metric components with respect to the (ro,Ti) coordinates are nearly constant, 
and thus this scaling still roughly holds. 

From the scaling of the proper volume we can obtain the scalings of the other quantities of interest. Inserting the 
dependence ( 2.12D of the proper volume per template on the minimal match, we obtain 




'^ — 1 ^— 8/3 ^ — 8/3 22 ^ 

Taking Eq. iH), neglecting the weak logarithmic dependence, and noting that fu is proportional to /o, we find that 

(3.13) 

Multiplying the number of templates Af by the length of each template F, the storage requirement scales as 

5 ^ (1 - (3.14) 


This last scaling with mmin presents a significant obstacle to efficiently searching for binaries composed of low-mass 
objects, such as MACHOs if they are low-mass (~ .5Mq) black holes [Q. Searches for low-mass objects will likely 
need to generate templates as needed rather than store them, incurring additional CPU costs not addressed in this 
paper. 

Note that the above scaling laws implicitly assume that /q or is varied while holding the overall shape of the noise 
spectrum fixed. The noise spectrum of a real interferometer, being composed of many independent noise sources, is 
unlikely to change in such a manner except for fairly small changes in /g or fg. However, these scaling laws give one 
a rough feel for (i) how changes in a single interferometer can affect data analysis requirements, and (ii) why different 
interferometers have drastically different requirements. 


IV. A TEMPLATE PLACEMENT ALGORITHM 

In this Section we deal with the actual placement of the templates in the parameter space. The parameter space 
is shown in Fig. ^ja in terms of the total mass M and symmetric mass ratio rj and in Fig. |b in terms of chirp times 
tq and Ti, for searches in which the maximum total mass M^ax of the system is 100 Mq and the lower limit Wmin 
on the mass of each component star is 0.2, 0.5 or I.OMq. The bottom line in Fig. ^ corresponds to binaries of equal 
mass {t] = 1/4) with the rightmost point corresponding to lowest mass binaries and the leftmost to greatest mass 
binaries of our search. There are no binaries in the region below this line as the parameter rj exceeds 1 /4 there, which 
is unphysical. The top and the left lines are determined by m-min and Mmax respectively. Given the minimum mass of 
the component stars and the maximum total mass, the parameter space of binaries is completely fixed. The volume 
of the parameter space (and the corresponding number of templates required) is a sharp function of the lower cutoff 
in the masses of the component stars and increases, as we have seen, as 

The shape of the parameter space is rather complicated and attention needs to be paid in the placement of templates 
so that the inevitable spill over (see below) is minimal. Our implementation of the filter placement is motivated by 
the following astrophysical consideration: The observed neutron stars are all of roughly equal mass Q. It is therefore 
to be expected that many inspiral signals will come from equal mass binaries. Consequently, we optimize the filter 
placement for equal mass binaries. This is achieved by beginning our template placement along the rj = 1/4 line. 
The span of each template is taken to be the largest rectangle (in a coordinate system in which the metric is locally 
diagonal) that can be inscribed inside the minimal match ellipse (we take MM = 97%). Note that in Fig. ^ the 
spans do not appear rectangular because they are sheared by transforming from coordinates in which the metric is 
locally diagonal (see below) to the (to,ti) coordinates. We begin with the leftmost point on the bottom edge of the 
parameter space of Fig. ^b, corresponding to the most massive binary of our search with the shortest chirp time. Such 
a system will, of course, consist of equal mass bodies and therefore our starting point is on the ?; = 1/4 curve. The 
next template is placed on the rj = 1/4 curve at that location where the left edge of its rectangle touches the right 
edge of the previous template’s rectangle. In a sense this is a straightforward generalization of placement of filters 
discussed in Ref. along grid lines that are not necessarily straight. 

This optimal translation of templates is most easily done in a coordinate system in which the metric is locally 
diagonal. Let (xo,xi) denote such a coordinate system (found by diagonalizing the two-dimensional matrix gij) and 
let f{xo,xi) = 0 denote a curve in the two-dimensional xq-xi plane along which templates are to be placed. For 
instance, 77 = 1/4 in Fig. || is one such curve. A convenient point to begin is the point at one end of the 
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FIG. 2. The parameter space of search in terms of (a) standard binary mass parameters and (b) chirp times, for different 


values of the lower mass limits. 











FIG. 3. Illustration of optimal translation of a template along an arbitrary curve. 


curve. In Fig. ^ we have sketched an arbitrary curve together with the first template and its span. The span of a 
template, for minimal matches close to I, is an ellipse. However, its effective (non-overlapping) span, when setting 
up a lattice of templates, is only an inscribed polygon such as a rectangle or an irregular hexagon. In the following 
discussion for simplicity we consider the span to be a rectangle and hence we will be setting up a rectangular lattice. 
By choosing a hexagonal lattice the number of templates can be reduced by about 20%, but the reduction is less when 
the curve along which templates need to be placed is parallel to neither Xg nor Xi axis. 

Given the ‘local’ distance (dxQ^\ dx[^^) between templates, we can get two points on the curve f(xo, Xi) = 0 which 
are simultaneous solutions of 


{/(a^o,a;i) = 0 , x, 


- 7.(1) 




(4.1) 


for n = 0 and n = 1 respectively. In order to cover the parameter space without leaving any ‘holes’ it is obvious 
that the next template should be placed at the point that is nearer to the first template. This is how one obtains the 
nearest neighbour of a template. 

Returning to our problem of placing templates in the parameter space of chirp times, we move along the 77 = 1/4 
line till the binary of longest chirp time is reached. The set of templates chosen on the 77 = 1/4 curve form the base, 
on top of which we construct layers to fill the parameter space. The next row of templates is set on top of the first 
so that the region of interest is completely covered by the spans—which means that the first few templates will be 
located outside that region, and that the total number of templates must be modestly greater than naively estimated 
by Eq. ( p.lO ) (see Fig. ^). 


V. DISCUSSION 

In this paper we have discussed the problem of optimally placing templates in a binary inspiral signal search. 
For our templates we have used the restricted second post-Newtonian waveform. We have made estimates of the 
computational costs of an on-line search of the inspiral waveform for all the large- and mid-scale interferometeric 
detectors now under construction. We have addressed several important issues: (i) the density of templates in the 
parameter space, (ii) the set of parameters most suitable for easy placement of templates, and (hi) the number of 
templates and computational resources needed to analyze the data on-line. 

These estimates should serve as a baseline for explorations of other data analysis strategies, and some of the 
techniques here could be incorporated into other strategies. For example, it is now recognized that a substantial 
reduction in computational cost can be achieved by carrying out a hiearchical search | p5|p^ . For a two-step hierarchical 
search strategy, in the first step a sparsely filled family of templates is used, with a threshold lower than what is 
acceptable based on the expected number of false alarms. Those events which cross this trigger threshold are further 
examined with a finer grid of templates chosen around the template that triggered the event. In such a hierarchical 
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Templates for M in the range [2,100] 



“o (s) 


FIG. 4. Choice of templates in the space of chirp times for initial LIGO interferometer with an upper frequency cutoff of 2 
kHz and for a search of binaries with masses of component stars larger than IMq and total mass no greater than IOOMq . The 
parameters characterizing the lattice at the location {M = lAMQ,r] = 1/4) are dxo = 46.5 ms, dx\ = 4.45 ms and 6 = 0.35 rad, 
and vary slightly from these values at other locations. 


search, templates chosen in the first step will essentially be the same for each data segment. However, templates 
in the second step need to be changed from one data segment to the next, depending on which templates from the 
coarse grid family produces an ‘event’. It is in the case of a coarse grid that our analytical algorithm for template 
placement and analytical estimates of computational requirements fail and must be replaced by numerical methods 
that are computationally expensive. But generating filters corresponding to the first step of the hierarchical search 
is more or less a one-time job. A finer grid is to be chosen quite frequently (essentially each time a possible event 
is selected in the first step) and in this case, fortunately, the analytical techniques discussed in this paper are quite 
accurate and one does not have to follow the time consuming numerical placement of templates either for estimating 
computational costs or for actually performing the search. 

There are several important problems we have not addressed in this paper which could be the topics of future 
work. We have not treated the problem of searching the corners of parameter space where the precessing binaries 
live. Although requiring some effort to seek out, these systems could prove quite informative and astrophysically 
interesting. The problem of searching for precessing binaries has been addressed only in a very exploratory way |lq ] 
but could benefit from further analysis using the techniques of this paper. Also, we have used a crude relation between 
the minimal match and fraction of event rate lost. This could be improved by a statistical analysis such as begun by 
Mohanty p^ . Now that the “P-approximants” [ p^ have proven a promising way of building templates, it is important 
to examine the computational costs of using them to conduct a search. 
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